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Abstract 

We express Maxwell's equations as a single equation, first using the divergence of a special 
type of matrix field to obtain the four current, and then the divergence of a special matrix to 
obtain the Electromagnetic field. These two equations give rise to a remarkable dual set of 
equations in which the operators become the matrices and the vectors become the fields. The 
decoupling of the equations into the wave equation is very simple and natural. The divergence of 
the stress energy tensor gives the Lorentz Law in a very natural way. We compare this approach 
to the related descriptions of Maxwell's equations by biquaternions and Clifford algebras. 



1 Introduction 

Maxwell's equations have been expressed in many forms in the century and a half since their dis- 
covery. The original equations were 16 in number. The vector forms, written below, consist of 
4 equations. The differential form versions consists of two equations; see [Misner, Thorne and 
Wheeler(1973); see equations 4.10, 4.11]. See also [Steven Parrott(1987) page 98 -100 ] The ap- 
plication of quaternions, and their complexification, the biquaternions, results in a version of one 
equation. William E. Baylis (1999) equation 3.8, is an example. 

In this work, we obtain one Maxwell equation, (10), representing the electromagnetic field as a 
matrix and the divergence as a vector multiplying the field matrix. But we also obtain a remarkable 
dual formulation of Maxwell's equation, (15), wherein the operator is now the matrix and the field 
is now the vector. The relation of the four vector potential to the electromagnetic field has the 
same sort of duality; see equations (13) and (14). 

These dual pairs of equations are proved equivalent by expanding the matrix multplication and 
checking the equality. Indeed, they were discovered this way. However, it is possible to ask if there 
is an explanation for this algebraic miracle. There is. It follows from the commutivity of certain 
types of matrices. This commutivity also implies the famous result that the divergence of the 
stress-energy tensor T is the electrom magnetic field F applied to the the current-density vector. 

The commutivity itself, again discovered and proved by a brute force calculation, has an explanation 
arising from two natural representations of the biquaternions; namely the left and the right regular 
representations. 



Why do matrices produce an interesting description of Maxwell's equations, when tensors are so 
much more flexible? It is because matrices are representations of linear transformations for a given 
choice of a basis. The basis is useful for calculation, but reformulating concepts and definitions 
in terms of the appropriate morphisms (in this case, linear transformations) almost always pays a 
dividend, as we topologists have discovered during the last century. 

This paper arose out of an email from me to Vladimir Onoochin concerning questions about 
Maxwell's equations. I thank him for a very interesting correspondence. 



2 Maxwell's Equations 

We are not really concerned with physical units for our paper, however what we have written is 
compatible with natural Heavyside-Lorentz units where the speed of light c = 1 and the electric 
permittivity eo = 1 . See Baylis(1999), section 1.1. 

We will say that two vector fields E and B satisfies Maxwell's equations if 



Here =: means the right symbol is defined by the left side of the equation. Now we follow [Got- 
tlieb (1998)] and [Gottlieb (2001)] and recall the notation for Lorentz transformations. Let M be 
Minkowski space with inner product ( , ) of the form — h ++. Let eo, e±, e 2 , e 3 be an orthonormal 
basis with eo a time-like vector. A linear operator F : M — > M which is skew symmetric with 
respect to the inner product ( , ) has a matrix representation, depending on the orthonormal basis, 
of the form 



V x E + d t B 

V x B — <9 t E 

V-E 



V B 
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J 
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(3) 
(4) 




(5) 



where xB is a 3 x 3 matrix such that (xB)v = v x B, the cross product of v with B. That is 



/ B 3 -B 2 \ 

-B 3 B 1 
V B 2 -B 1 / 



where B is given by {B\, B 2 , B 3 ). 
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The dual F* of F is given by 



F* : = ( °. ~ B ! ) (6) 
\ -B xE y ' 



We complexify i 7 by cF := F — iF* . Its matrix representations is 

cF = \ ° r , A - ] where 1 = £ + iB (7) 
^ A x(-tA) j 

Remark 1 These complexified operators satisfy some remarkable Properties: 



a) cF\cF2 + CF2CF1 = 2{A\,A2)I where { , } denotes the complexification of the usual inner 

product of M 3 . Note, { , } is not the Hermitian form, that is our inner product satisfies 
i(v,w) = (iv,w) = (v,iw). 

b) The same property holds for the complex conjugates cF\ and C.F2 . The two types of matrices 

commute. That is cF\cF2 = CF2CF1 

c) cFcF =: 2Tp where Tp is proportional to the stress-energy tensor of electromagnetic fields. 



Now e F = I + F + F 2 /2\ + i ?3 /3! + . . . : M — > M is a proper Lorentz transformation. It satisfies 
e F = e cF l 2 e cF l 2 . The algebraic properties of cF give rise to a simple expression for its exponential: 

cF u/\ \t , sinh(A c ,F) 
e = cos1i(A c _fK H -cF 

where X c f is an eigenvalue for cF. There are only, at most, two values for the eigenvalues of cF, 
namely X c f and —X c f- 



Now let us see how to describe Maxwell's field equations by means of these matrices. Now consider 
the 4x4 matrix F defined by equation ©. We multiply it by the 4x1 vector of operators (— dt, V)* 
to get 

F ( "v) - G p ) < 8 > 



and 



-<v'H:) 



See equation (23) for more details on this notation. 
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Equations |H] and hold if and only if Maxwell's equations (1-4) hold. Using the definition of cF 
in equation (7) we get 



-I v ) = <-) 

is true if and only if E and B satisfy Maxwell's equations. 

In fact, the version of Maxwell's equations involving the four vector fields E, B, D, H is equiv- 
alent to the matrix equation 



D + iB \(-9t\ = (p\ (n) 

D + iB x(-i(E + iH)) J V V J \jj 1 j 



Remark 2 Note the forms of the complex matrices in equation (7) and in equation (11): respec- 
tively 



( A x(±L4) ) and ( A x(C) ) 

Now within the class of all matrices of the second form, only the first form has the property that 
its square is equal to a multiple of the identity. 



Now let us consider a four vector field f j? J . We define an associated E and B by the equations 

V x A =: B and - <9 t A — V<p =: E (12) 



Then the E field and the B field defined above satisfy Maxwell's equations (1 - 4). We can describe 
the four vector field by means of a similar matrix equation. Let / denote the 4x4 identity matrix. 



Now a remarkable duality holds. The following equation is also equivalent to the above equation. 
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Similarly, equation (10) has a dual equation which holds if and only if Maxwell's equations are 
satisfied, 



«-i:,r v ))u>(;) w> 



These equations give rise to an interesting question: 



Remark 3 Can we give an explanation of the remarkable dualities between equations (10) and 
(15) and between equations (13) and (14) ? Yes! It is based on the fact that matrices of the form 
cF commute with the matrices of the form cF. Hence the matrix on the left hand side of (13) 
commutes with the matrix on the left hand side of of (14)- Thus the first columns of two products 
of the commuting matrix products must be equal. Since the vectors on the left hand side of (13) 
and 14 ) are the first columns of each of the two product matrices that commute, it follows that their 
products are equal. Exactly the same argument holds for equations (10) and (15). 



Now from part a) of Remark 1, we know that 



V* 

V ixV 



V* 

V ixV 



2 

2 



= V 2 / (16) 



Consistent with our notation we define 



=: cV (17) 



so that 

(dti - cV){d t i + cV) = (d t i + cv)(d t i - cV) = (d 2 - v 2 )/ (is) 

Now apply (dtl + cV) to equation (15) and obtain the wave equation and obtain the wave equation 

(d 2 - v 2 ) (_ E °_ . B ) = Vp + d t J - iV x J (19) 

implying the conservation of charge equation dtp + V • J = 

On the other hand, if we apply [dti — cV) to equation (14) we obtain, thanks to equation (15), 

(<9 2 - V 2 ) (j^j = (^j + vector depending on d t <p + V ■ A (20) 

So if the covariant gauge condition is chosen, dt<p + V • A = 0, we have the wave equation 

CO -ft) < 2i > 
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3 The electromagnetic stress-energy tensor 



Not only is equation (10) equivalent to Maxwell's equations, but its complex conjugate below is 
also equivalent to Maxwell's equations, since the current is real. 



Thus we have equations (10), its complex conjugate (22), and its dual (15), and the complex 
conjugate of (15) all being equivalent to the vector form of Maxwell's equations (1-4). Thus we 
not only reduce the number of Maxwell's equations from 4 to 1, we obtain 4 equivalent equations. 
Thus, by a talmudic argument, we can say we have reduced the four Maxwell's equation to 1/4 
of an equation. These different forms of the equation interact with each other to produce new 
derivations of important results. 

We have used the following notation of the divergence of a matrix field. What we mean by the 
notation 



is that the column vector of operators multiplies into the matrix A of functions and then the 
operators are applied to the functions they are next to. In index notation we obtain a vector whose 
i-th row is ciijdj := dj(aij). 

Another way to achieve the same result is to take the differential of the matrix, dA. Here we obtain 
a matrix of differential one-forms whose (i, j)-th element is 



Then we employ the differential geometry convention that dxi{dj) = 5ij. Thus da,ij(dj) = dj(aij). 
Thus our definitions leads us to 




(22) 




(23) 



dciij = dk{aij)dxk 



(24) 




(25) 



Now the Leibniz rule gives us 



d{AB) = (dA)B + A(dB) 



(26) 
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This equation helps us to study the divergence of a matrix product. In particular, if dA and B 
commute, then the divergence of the product is B times the divergence of A plus A times the 
divergence of B. 



Theorem 1 Let Tp be the electromagentic stress energy tensor of the electromagnetic matrix field 
F. Then 

-<9A 



Proof: 



The first equation below follows from (25). The second equation follows from Remark 1 c). 
The third equation follows from (26). The fourth equation follows from Remark 1 b), which 
is the commutivity of matrices of the form cF with matrices of the form cG, combined with the 
observation that d(cF) has the same form as cF and similarly d(cF) has the same form as cF. 
The fifth equation follows by linearity . The sixth equation follows from the form of Maxwell's 
equations found in equation (10) and its complex conjugate equation (22). Finally the seventh 
equation follows from the definition of cF and the fact that F is its real part. 

T F = d{T F ) ■ = \d(cFcF) ■ = l -(d(cF)cF + cFd(cF)) • 

\{cFd{cF) + cFd(cF)) ■ = \cFd{cF) ■ + \cFd(cF) • 

H?)+H;)-'(; 

□ 



4 Biquaternions 

We choose an orthogonal coordinate system (t,x,y,z) for Minkowski space R 3 ' 1 . Let be the 
corresponding unit vectors (with respect to the Minkowski metric). We express e± = (1, 0, 0)', e*2 = 
(0, 1, 0)*, e*3 = (0, 0, 1)*. So our standard choice of basis is given by e$ = (0, e^)* for i ^ 0. 

Let cEi be the matrix below where A = e*j where i = 1, 2, 3. 



cE i =\ r , A where A = E + i5 = ei (28) 
V ^ x(-tA) / 



Remark 4 aj T/ie set of sixteen matrices 

I, cEi, cEi, cEicEi 
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where i = 1,2, 3, forms a basis for M^(C), the vector space of 4 x 4 complex matrices. 

b) The square of each matrix in the basis is I. 

c) Each matrix is Hermitian, and all the matrices, except for the identity I, have zero trace. 

d) cE\cE2 = icE% and cEicE 2 = —icE 3 . 

e) cEi and cEj anti commute when i ^ j; and cEi and cEj anti commute when i ^ j. Also cEi 
and cEj commute for all i and j. 

See Gottlieb(2001), Theorem 3.3. 

Now from the above theorem, we see that the 4 elements {/, cEi} form a basis for the biquaternions. 
Baylis(1999) chooses to denote his basis for the biquaternions, which he views as the real Clifford 
algebra on three generators, with a basis ei where i = 0, 1, 2, 3 where eo is the multiplicative identity 
and e? = 1, and the ei anti commute for i ^ 0. Baylis chooses the orientation eie2 = ie3. Thus 
our representation of Baylis's biquaternions is given by eo i— ► / and e» i— > oE^. 

Let us denote the biquaternions given above by P. The alternative choice where the biquaternions 
satisfy the relation eie2 = — ie3 will be denoted P and is represented by e, i— > cEi. The symbol P 
stands for paravectors, which is Baylis' name for the space of biquaternions. The name paravectors 
stems from Baylis' underlying Clifford algebra approach. 

Given a paravector A £ P, we define left multiplication by A as 

L A : P -» P : X ^ AX 

and right multiplication by A as 

R A ■ P -> P : X ^ XA 

Lemma 2 Both left multiplication and right multiplication by A are linear transformations, and 
they commute: that is LaRb = RbLa 

Now given the basis {ei}, the linear transformations can be represented by matrices muliplying the 
coordinate vectors from the left. 

Theorem 3 L A is represented by the matrix al + cF corresponding to A. Also, Ra is represented 
by the transposed matrix (al + cF)*. If A G P, then A is represented by al + cF and so Ra is 
represented by the transposed {al + cF) t . 

proof: Consider L\, by which we mean left multiplication by ei. So the basis elements are trans- 
formed by L\\ eo >— > eieo = ei, eo i— > eiei = eo, eo >— > eie2 = ie3, eo i— > eie3 = — ie2- This 
corresponds to the matrix 
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/ 1 \ 

10 

-i 

\ i / 



This matrix is cE\. Right multiplication by ei gives us the matrix cE\. Note that this is the 
transpose of cE\. In the same way we see that left multiplication by ei gives rise to the matrix cE{ 
and right multiplication by ei gives rise to cEj = cEi. 

Now every element A € P is a unique linear combination a^ei, so the representation of A is the 
same linear combination of the cEi : Namily, aicEi. Thus La is represented by ciicEi. Also Ra is 
represented by (ajcE'j)* = a,jcEi □ 

Baylis(1999) in his textbook, describes an algebraic system in which the paravectors are uniquely 
written as linear combinations a^ei. The choice of the basis implicitly defines an isomorphism 



Whereas if v is a nonnull vector in complex Minkowski space M the evaluation map is an isomor- 
phism 



proof: See Gottlieb(1999), Theorem 6.9. 

Note that if v = eo, then <E> = eo . This follows since aicEi i— > aicEieo = a%ei = (ao, ai, ct2, azf . 

Now Baylis(1999) exposes Electromagnetism in the language of Clifford algebra on three generators, 
which is isomorphic to the biquaternions. We therefore want to adopt his notation as much as 
possible, while thinking of our matrices replacing his abstract symbols. Our goal is not to advance 
a new method of calculation so much as to use the matrices to understand the underlying geometry 
which arises out of an electromagnetic field. Our point of departure, explained in Gottlieb(1999), 
is to consider the linear transformation F : M — > M which arises in the Lorentz law. Choosing 
an orthonormal basis {ej} with respect to the Minkowski metric <, >, we consider those linear 
transformations F : M — ► M which are anti symmetric with respect to the Minkowski metric. 



$ : P — > C 4 : X = a^ei i-> (a , ai, a 2 , a 3 ) 



t 



(29) 



0^ : P -> C 4 : X = aicEi ^ Xv 



(30) 



That is 



F = -F 



(31) 



where F is defined by the equation 



< Fv, w >=< v, Fw > 



(32) 



That is, F is the adjoint of F with respect to the Minkowski metric. 
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The matrix form which these skew symmetric operators take is that of the matrix in equation 
(5). Now our point of view is the following. The Minkowski metric is taken to be the primary 
object, based on the work of Robb(1936), who produced a protocol as to how the metric could be 
measured by means of light rays. Also we take linear transformations as primary, since they are 
the morphisms of the category of vector spaces and linear transformations. This point of view is 
suggested by the successes of algebraic topology. We do not hesitate to employ other inner products 
on M or the convenience of tensors, but we will never change the underlying sign of the metric to 
match the sign suggested by the wave equation. 

For example, the Minkowski metric <, > is taken to be defined independently of a choice of basis. 
Having chosen an orthonormal basis, we of course represent F by a matrix of type (5). We can also 
define based on this choice of a basis, the Euclidean metric <, >c and the Hermitian metric <, 
Now the adjoints of F are defined by < F t v,w >c*=< v,Fw >c and < F^v,w >h=< v,Fw >h- 
On the matrix level, t is the transpose and f is the complex conjugate transpose. The adjoint 
{al + cF) = al - cF. 

Now Clifford algebras have certain involutions which are defined ad hoc and are very useful in 
Clifford algebras. These are Clifford conjugation (or spatial reversal) denoted by p i— > p and 
hermitian conjugation (or reversion) denoted by p i— > p' 

Now starting from a real linear transformation, the definition of cF is seems just a clever trick that 
is very useful. If you take a point of view that the biquaternions are the key concept, then the 
complexification trick is explained. Similarly, the Clifford conjugation and hermitian conjugation 
are tricks from the point of view of Clifford algebra, but from the matrix point of view they 
correspond to the Minkowski adjoint and the Hermitian adjoint respectively. 

Another important involution is complex conjugation denoted by p i— > p c The bar is reserved for 
Clifford conjugation, and the * is reserved for the Hodge dual , see equation (6). Now the space 
of 4 x 4 matrices is a tensor product of P and P. Now transpose t : P — > P and complex 
conjugation c : P — > P. The composition of complex conjugation and transpose is the Hermitian 
adjoint f : P — > P. Both the Minkowski adjoint and the Hermitian adjoints reverse the orders of 
multiplication because they are both adjoints. Thus AB = B A and (ABy = B^A^. Complex 
conjugation is important on the underlying vector space C 4 where it is given in our new notation 
by (ajej) c = afej. Using the isomorphism <E>, see (29), it can be grafted onto the biquaternions as 
(ajei) c = a?ei. This definition is not used in Baylis(1999) probably because it has no simple relation 
to the ring structure of the biquaternions. However, the natural extension of complex conjugate to 
matrices, given by (aij) c = (a?-)> is important on the matrix level. For our matrix description of 
biquaternions, we get 

(al + cFf = a 5 1 + cF (33) 

If we apply transpose to this equation we get 

(al + cF) 1 " = a? I + cF ] (34) 

The problem for a strictly biquaternion approach is that the right side of (33) is no longer in the 
the biquaternions P, and although c does preserve order of products, it also is not linear over the 
complex numbers. 
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The composition of Hermitian conjugation with Clifford conjugation leads to an product order 
preserving automorphism + : P — > P : (al + cF) i— » (al + cF) + = a c I — cF^ . This automorphism 
is called the grade automorphism. 

Our Hodge dual, see (6), eventually differs with that of Baylis. We have two choice for the Hodge 
dual, either the matrix in (6) or its negative. With our definition cF := F — iF* we agree with 
Baylis' choice of orientation e\e2 = iej, via the left regular representation. Then extending the 
definition of dual to cF we use cF* = (F — iF*)* := F* - iF** = F* +iF = i(F - iF*) = icF. So 
in P it looks like the definition should be multiplication by i. If we used the alternative definition 
of (cF)* = —icF* this gives multiplication by —i. But P agrees with the orientation e\e2e^ = i. 
Baylis chooses to define the Hodge dual by p* := pe^ei, which is effectively multiplying by —i. 



5 Biparavectors 

Let cEq := I, and let cEq := I Then {cEi <g> cEj} forms a vectorspace basis for P <g> P. Similarly, 
{cEi <g> cEj} forms a vectorspace basis for P <g) P. Now we will call any element in P <g> P a 
biparavector. So a biparavector is a sum of 16 terms give by aijcEi ®cEj. Now a bivector gives rise 
to a linear transformation T : P — > P : X i— > aijcEiXcEj In fact, every linear transformation can 
be represented this way as a biparavector. Now by Theorem 3 the linear transformation given by 
the biparavector {a^cEi ® cEj} is represented by the matrix {aijcEicEj} 

Theorem 4 Let A = acijEi <g> cEj be a biparavector representing a linear transformation M, of 
C 4 , repersented by the matrix again called M. 

a) Then M = aijcEicEj where aij = jtrace(McEicEj). 

b) Let x £ C 4 correspond to X € P by Xeo = x. Then Mx = ^(J^eo := aijcEiXcEjeo 

proof: follows from remark 4. The trace from part a) follows from Remark 4, c). The trace of 
the product of two matrices tr(AB) = tr(BA) gives rise to an innerproduct on the vector space of 
matrices. The the above cited remark shows that the basis in question is orthnormal with respect 
to the inner product given by \trace. 

Corollary 5 a) The stress-energy tensor Tp = \cFcF , so the equivalent biparavector is cF ® cF^ 
b) A real proper Lorentz transformation e F = e^ cF e^ cF So the equivalent biparavector is e^ cF ®e^ cF<i 

proof: The key point is that cF* = cF^ 
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